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摘 要 : 为 探究 黄河 上 游 径流 量 与 降水 量变 化 特征 ,基于 唐 乃 辫 、 下 河沿 、 头 道 描 3 个 水 文 站 近 60 a 
的 降水 量 、 径 流量 资料 ,采用 线性 倾向 趋势 检验 、Mann-Kendall 检 验 法 (M-K 趋 势 检 验 )、Spearman 秩 
次 相关 检验 、M-K 突变 检验 、Pettitt 非 参数 检验 、 有 序 聚 类 分 析 、 累 积 距 平 法 、 双 辕 积 曲线 等 分 析 方 


法 对 比 研究 了 黄河 上 游 及 上 游 不 同 子 区 域 降水 量 和 径流 量 的 变化 特征 ,并 讨论 了 径流 量 对 降水 量 


的 响应 关系 。 结 果 表 明 : 黄 河上 游 降水 量 呈 不 显著 增加 趋势 ,在 2003 年 发 生 突变 ,突变 前 后 变化 率 
为 4.67%; 径 流量 呈 显 著 减 少 趋势 ,突变 年 份 为 1986 年 ,突变 前 后 变化 率 为 33.34%。3 个 子 区 域 年 
降水 量变 化 趋势 分 别 呈 显著 增加 、 不 显著 增加 和 显著 减少 趋势 ,径流 量 均 呈 减少 趋势 。 以 唐 万 记 
以 上 的 工区 为 基准 区 时 ,上 唐 乃 亥 一 下 河沿 所 在 的 下 区 降水 量 因素 对 径流 量 的 影响 达到 25.08% , 非 


降水 量 因素 为 74.92% , 而 下 河沿 一 头 道 拟 的 亚 区 ,降水 量 对 径流 量 的 影响 为 32.14% , 非 降 水 量 因 


素 为 67.86%。 研 究 结 果 对 黄河 流域 水 资源 综合 管理 与 科学 配置 具有 参考 意义 。 
X 键 词 : 黄河 上 游 及 分 区 ; 径流 量 ; 降水 量 ; 变化 特征 ; 响应 关系 


文章 编号 : 


降水 量 与 径流 量 作为 自然 界 水 循环 的 2 个 重要 
环节 对 区 域 水 资源 产生 重大 影响 ,降水 量 是 区 域 水 
资源 量 时 空 分 布 特征 最 重要 的 因素 ,也 是 区 域 自然 
地 理 特征 重要 表征 要 素 与 关键 环 市 ;径流 量 是 水 资 
源 的 重要 组 成 部 分 和 主要 存在 形式 ,是 地 表 水 循环 
的 重要 环节 和 水 量 平衡 的 基本 要 素 ,更 是 社会 经 济 
发 展 用 水 的 最 主要 来 源 。 然 而 ,径流 量 的 变化 受降 
水 量变 化 的 影响 ” ,因此 ,研究 区 域 降水 量 径流 量 
变化 规律 与 特征 以 及 二 者 之 间 的 关系 和 影响 因素 ， 
掌握 区 域 水 文 水 资源 演变 规律 ,可 为 区 域 水 资源 综 
合 管理 与 科学 配置 提供 重要 参考 和 借鉴 依据 。 黄 
河 是 中 国 第 二 长 河 ,是 中 国 西北 、 华 北 地 区 的 重要 
水 源 ” 。 黄 河流 域 西 起 巴 颜 喀 拉 山 , 东 临 渤海 , 南 至 
秦岭 , 北 抵 阴 山 ,是 横 跨 青藏 高 原 .内 蒙古 高 原 、 黄 
土 高 原 和 黄 淮海 平原 的 生态 廊 道 ,承载 着 全 国 1596 
的 耕地 面积 和 12% 的 人 口 扩 ,沿线 集聚 有 关中 平原 、 
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中 原 .山东 半岛 3 个 区 域 级 城市 群 和 兰 西 .宁夏 沿 
黄 . 呼 包 哪 榆 、 晋 中 4 个 地 区 性 城市 群 守 , 在 中 国 经 
济 社 会 发 展 和 生态 安全 中 占有 十 分 重要 的 地 位 。 
但 黄河 流域 受 地 理 位 置 与 东亚 季风 气候 的 影响 , 流 
域 大 部 分 属于 干旱 、 半 干旱 地 区 ,多 年 平均 降水 量 
在 200~600 mm 之 间 , 水 资源 自然 课 赋 严重 不 足 “， 
人 均 水 资源 占有 量 仅 为 全 国平 均 水 平 的 27% 。 
2016 年 ,流域 缺 水 8.60x10" m: , 缺 水 率 高 达 17% (全 
HKKK 2.2%)" ,供需 矛盾 十 分 突出 。 

因此 ,众多 学 者 聚焦 影响 黄河 流域 水 资源 关键 
水 文 要 素 ( 径 流量 与 降水 量 ) 及 影响 因素 展开 了 深 
入 而 广泛 的 研究 ,如 保 广 裕 等 利用 气象 观测 站 逐 
日 地 面 降 水 量 和 水 文 站 径流 量 数据 ,分 析 了 黄河 上 
游 河 源 区 不 同 量 级 降水 量 对 径流 量变 化 的 影响 ; 魏 
伊 宁 等 ”研究 了 兰州 以 上 区 工人 径流 量 丰 枯 的 空间 分 
布 特征 及 其 影响 因子 ; 宁 怡 楠 等 分析 黄河 中 游 河 
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成 艺 等 : 黄河 上 游 近 60 a 径流 量 与 降水 量变 化 特征 研究 


龙 区 间 4 个 典型 流域 径流 量变 化 趋势 ,阐述 了 气候 变 
化 和 人 类 活动 对 流域 径流 量变 化 的 作用 ; 潘 彬 等 
对 黄河 下 游 1962 一 2012 年 径流 量 和 降水 量变 化 特 
征 及 规律 进行 了 深入 探讨 ; 王 未 等 ”利用 多 维度 径 
流量 情势 指标 ,分 析 了 黄河 流域 各 区 域 的 径流 量 过 
程 整体 变化 特征 并 评估 了 其 改变 度 。 尽 管 国内 外 
学 者 对 黄河 流域 降水 量 径流 量变 化 规律 和 影响 
素 采 用 不 同方 法 .不 同 站 点 数据 及 不 同时 间 序 列 多 
有 探究 ,但 对 降水 量 径流 量 时 空 尺度 变化 特征 研 
究 的 成 果 存在 一 定 差异 后 ,针对 流域 上 游 降水 量 、 
径流 量 的 研究 大 多 集中 在 以 唐 乃 北 水 文 站 以 上 的 
河源 区 ,涉及 整个 上 游 区 域 (河源 一 头 道 拐 ) 的 研究 
较为 单薄 ,上 且 绝 大 多 数 人 研究 都 集中 在 时 间 序 列 突变 
前 后 二 者 的 变化 及 影响 因素 的 分 析 ,而 针对 上 游 分 
区 域 降水 量 、 径 流量 响应 的 研究 几乎 空白 。 因 此 ， 
本 文 根 据 黄河 上 游 地 区 地 域 特征 ,以 典型 水 文 站 为 
界 , 将 黄河 上 游 分 为 3 个 研究 区 域 ,分 区 域 研 究 黄 河 
上 游 降水 量 、 径 流量 变化 趋势 及 二 者 之 间 的 响应 ， 
为 黄河 流域 生态 保护 与 高 质量 发 展 提供 科技 支撑 ， 
夯实 黄河 全 流域 高 质量 发 展 根基 。 


1 研究 区 概况 


根据 黄河 流域 形成 发 育 的 地 理 .地质 条 件 及 水 
文 情 况 ,将 源头 至 内 蒙古 自治 区 托 克 托 县 河口 镇 区 
域 称 为 黄河 上 游 ,该 区 域 ( 含 内 流 区 ) 位 于 青藏 高 
JE .内 蒙古 高 原 和 黄土 高 原 三 大 高 原 交接 地 带 , 流 
域 面积 3.86x10’ km , 占 流域 总 面积 的 53.80% ,径流 
量 占 全 河 的 60% 以 上 ,如 单纯 选择 整个 黄河 上 游 
进行 研究 并 不 能 很 好 地 反映 区 域 降水 量 径流 量变 
化 的 时 空 规律 。 因 此 ,基于 水 资源 二 级 区 正式 图 
件 ,综合 流域 自然 环境 和 人 类 活动 对 降水 量 .径流 
量 影响 的 差异 ,以 流域 典型 水 文 站 断面 为 界 ,将 研 


头 道 拟 水 文 站 之 间 的 区 域 ,平均 海拔 1344.09 m, E 
年 平均 气温 7.00 % ,多 年 平均 降水 量 271.7 mm, 区 
域 地 势 平坦 ,分布 着 广 亡 的 耕地 资源 ,是 黄河 上 游 
主要 的 灌溉 农业 区 ,该 区 域 的 宁 蒙 河套 灌区 是 全 国 
特大 型 灌区 之 一 , 赖 于 黄河 水 的 灌溉 ,河套 平原 享 
有 “ 八 百 里 河套 米 粮 川 " 的 美誉 (图 1)。 


2 数据 与 方法 


2.1 数据 来 源 

本 研究 所 用 径流 量 数据 为 1956 一 2017 年 黄河 
上 游 唐 乃 北 、 下 河沿 及 头 道 抛 3 个 水 文 站 年 实测 径 
流量 数据 ,数据 来 源 于 国家 科技 基础 条 件 平台 国家 
地 球 系统 科学 数据 中 心 (http: /www.geodata.cn )。 研 
究 区 区 径流 量 数据 为 唐 力 北 水 文 站 实测 数据 ; 
IL HI 区 采用 区 域 径流 量变 化 量 , 分 别 为 相应 水 文 
站 ( 唐 力 辫 一 下 河沿 、 下 河沿 一 头 道 抛 ) 径 流量 差 值 
数据 ;上 游 径流 量 数 据 采 用 头 道 抛 水 文 站 实测 数 
据 。1956 一 2017 年 降水 量 数据 是 经 496 个 独立 气象 
观测 点 进行 验证 的 可 信和 数据 "1。 
2.2 研究 方法 

(1) 趋势 检验 

线性 倾向 估计 是 目前 趋势 分 析 中 最 直接 有 效 
地 方法 之 一 ,建立 分 析 要 素 % 与 时 序 i 之 间 的 一 元 函 
数 ,来 检验 时 间 序 列 的 趋势 性 "”。 滑 动 平均 是 趋势 
拟 合 技术 最 基础 的 方法 之 一 ,此 方法 在 降低 数据 序 
列 自由 度 的 同时 能 增强 序列 的 相关 性 ,可 采用 时 间 
序列 的 平滑 值 来 显示 变化 趋势 ,从 滑动 平均 曲线 图 
来 判别 其 变化 是 呈 上 升 或 下 降 趋势 。 

Mann-Kendall 检验 法 (简称 M-K 趋 势 检 验 ) 是 一 
种 非 参数 检验 方法 ,优点 为 样本 序列 无 需 服 从 特定 
分 布 , 且 不 易 受 奇异 值 干扰 ,用 于 水 文 要 素 时 间 序 
列 的 趋势 检验 。 此 方法 利用 统计 量 Z 值 的 正 负 判断 


究 区 划分 为 3 个 子 区 域 : E 区 为 河源 一 唐 乃 亥 水 文 站 之 
间 的 区 域 ,平均 海拔 4063.23 m, 多 年 平均 气温 -2.60%C， 
多 年 平均 降水 量 537.9 mm, 为 黄河 上 游 主 要 产 水 区 , H. 
位 于 连续 多 年 冻 土 区 ,多 年 冻 土 占 该 区 域 面积 的 
34%。 了 区 为 唐 力 辫 一 下 河沿 水 文 站 之 间 的 区 域 ， 
平均 海拔 2709.47 m, 多 年 平均 气温 3.08 %C ,多 年 平 
均 降 水 量 449.3 mm, 为 高 山 与 峡谷 串珠 状 连 接地 
带 , 并 有 多 条 河流 汇 信 ,水 能 资源 丰富 ,分布 着 诸如 
龙 羊 峡 、 刘 家 峡 等 大 型 水 利 枢纽 。 亚 区 为 下 河沿 一 


序列 趋势 , 当 Z 为 正 值 时 ,序列 呈 上 升 趋势 ,反之 为 
下 降 趋 势 ,Z 的 绝对 值 越 大 ,序列 趋势 变化 越 显 著 ， 
对 于 给 定 的 显著 性 水 平 w=0.05, "4 ZISAZ ERE, 1A OS 
序列 趋势 显著 1。 

Spearman 秩 次 相关 检验 是 一 种 非 参 检验 方法 ， 
该 方法 计算 简单 ,精确 性 高 ,根据 秩 次 相关 系数 (7.) 
判别 时 间 序 列 的 趋势 性 。 当 =*>0 时 ,时 间 序 列 呈 上 
升 趋势,r<0 时 , 则 旦 降低 趋势 ,并 用 1 检验 法 检验 趋 
势 的 显著 性 “”。 


202208.00065v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


Uu TEE! 


4 


青 
0_125 km — 刘 家 峡 水 电站 Y 


(a) 黄河 上 游 


图 例 


。 水 电站 -一 黄河 干流 
4 水 文 站 ”一 一 黄河 支流 
高 程 /m C 流域 边界 


" 6235.33 æm 湖泊 
983.89 


图 1 研究 区 概况 


Fig. 1 Overview of the study area 


(2) 突变 检验 

M-K 突 变 检验 是 在 M-K 趋 势 检 验 基 础 上 建立 的 
一 种 确定 时 间 序 列 突变 的 方法 ,根据 构造 的 秩序 列 
样本 计算 统计 量 和 如 ,根据 i 和 Ui 曲线 在 置信 
范围 内 的 交点 ,确定 时 间 序 列 的 突变 点 1。 

Pettitt 检验 为 非 参数 检验 方法 ,方法 结构 简洁 ， 
受 异常 值 影响 小 ”1。 统 计量 IU, rl 最 大 值 对 应 的 时 
间 点 4 即 为 时 间 序 列 的 突变 点 ,Ui 为 : 


N 
Lv = U, y+ S sen; 7x) (1) 
FE 


式 中 :UV 为 统计 量 ,其 本 质 为 前 1 个 样本 值 大 于 后 1 

个 样本 值 的 个 数 ;i 为 样本 的 排列 序号 ;NN 为 样本 总 

数 ; ox, vox; 为 相同 时 间 序 列 的 2 个 样本 值 。 对 于 突 
变 点 的 显著 性 根据 构造 的 统计 量 P 进 行 检验 : 

P=2exp[-6U, ,/(N? + N*)] (2) 

式 中 :PP 为 统计 量 ; 当 P 小 于 显著 水 平 a=0.05, 则 认 
为 检测 出 来 的 突变 点 在 意义 上 发 生 显著 变化 。 

有 序 聚 类 分 析 法 属于 非 参 数 检验 方法 ,其 实 

质 是 计算 最 优 分 割 点 ,使 同类 之 间 的 离 差 平方 和 较 

小 ,而 不 同类 之 间 的 离 差 平方 和 较 大 。 

S, (7) 7 min » 一 元 小 十 > (x; a (3) 

式 中 : S, 为 离 差 平 方 和 ;7 为 假设 序列 分 割 点 ; 元 为 

分 割 点 前 的 序列 均值 ; x_, 为 分 割 点 后 序列 均值 ;x 


为 序列 值 ;分 别 计算 + 前 后 序列 的 离 差 平方 和 
5S,(7) , 当 S, (c) 为 最 小 时 ,所 对 应 的 r 就 是 最 优 分 割 


点 , 即 最 可 能 的 突变 点 。 

累积 距 平 法 利用 累积 距 平 曲线 反映 变化 趋势 ， 
耕 累 积 距 平 曲线 显示 上 升 趋势 , 则 时 间 序 列 距 平 值 
增加 , 若 累 积 距 平 曲线 显示 下 降 趋 势 , 则 距 平 值 减 
少 , 同 时 可 利用 累积 最 大 绝对 值 判断 出 时 间 序 列 突 
变 发 生 的 大 致 时 间 。 

(3) 双 累 积 曲 线 法 

双 累 积 曲 线 法 是 水 文 气象 要 素 一 致 性 分 析 中 
最 常用 的 一 种 方法 ,同时 分 析 水 文 要 素 的 变化 规律 
及 趋势 。 如 果 未 受到 外 界 干扰 时 ,存在 以 下 关系 : 

YR-:&YP-b (4) 

式 中 : > RON BUE RCUR k NRA MP 为 累积 
降水 量 ;b 为 常数 。 将 受到 干扰 的 累积 降水 量 作 为 
变量 VP ,可 根据 式 (4) 计 算 受 到 影响 的 累积 径流 
ht > R ,模拟 和 实测 径流 量 深 的 差 值 则 为 径流 量 受 
非 降水 量 因素 而 产生 的 变化 量 。 


3 结果 与 分 析 


3.1 降水 量 与 径流 量变 化 特征 

3.1.1 黄河 上 游 降 水 量 与 径流 量变 化 趋势 ”对 黄河 
上 游 ( 头 道 拐 以 上 区 域 ) 的 年 尺度 降水 量 数据 及 头 
道 拐 水 文 站 实测 径流 量 数据 进行 分 析 ,根据 (图 
2a~p) 所 示 的 降水 量 、 径 流量 可 知 ,黄河 上 游 年 均 降 
水 量 约 为 419.9 mm,1956 一 1970 年 年 际 变化 波动 幅 
度 较 大 ,1970 年 以 后 较 小 ,最 大 降水 量 出 现在 1968 
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图 2 黄河 上 游 降水 量 与 径流 量变 化 特征 
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Fig.2 Characteristics of precipitation and runoff in the upstream of the Yellow River 


4E. , 73 559.3 mm, 最 小 降水 量 为 1966 年 的 318.3 mm, 
相差 约 241.0 mm, 极 值 比 为 1.76。 以 5 a 滑动 平均 过 
程 线 分 析 , 降 水 量 序列 在 1980 一 2002 年 左右 呈 减 少 
趋势 ,2002 年 以 后 呈 增 加 趋势 ,由 线性 倾向 趋势 可 
知 黄河 上 游 降水 量 以 1.80 mm * (102) H EREM 
增加 趋势 。 与 降水 量 相 比 ,径流 量 呈 明显 减少 趋 
势 ,减少 速 率 为 1.94x10" m.(10aj ,最 大 径流 量 出 现 
JE 1967 4E£(4.38x10" m) , 约 为 最 小 径流 量 出 现年 份 
1997 年 (1.12xl00m) 的 3.91 倍 。 从 Sa 滑动 平均 过 
程 线 看 出 ,径流 量 呈 现 多 段 上 升 一 下 降 一 上 升 变 化 
过 程 ,在 1986 一 2004 年 左右 呈现 显著 减少 趋势 。 

采用 M-K FI Spearman 秩 次 相关 检验 法 验证 
1956 一 2017 年 降水 量 和 径流 量变 化 趋势 及 其 在 置 
信 水 平 a=0.05 下 的 显著 性 ,结果 如 表 1 所 示 。 

根据 检验 结果 ,降水 量 的 M-K 检 验 统 计量 Z= 
0.78, H.17120.78«Z,571.96, Spearman 秩 次 相关 系数 
r=0.09>0, 且 统计 量 Im=0.78<Tw=2.01 ,说 明 两 种 检 
验方 法 均 得 出 黄河 上 游 近 60 a 来 降水 量 呈 不 显著 增 
加 趋势 ;径流 量 的 M-K 统 计量 Z=-3.82 ,上 且 |Z1=3.82> 
Zan=1.96 , Spearman 秩 次 相关 系数 =-0.51<0, 且 统 
计量 I7I=2.58>7w, 说 明 两 种 检验 方法 均 得 到 黄河 上 


游 近 60 a 来 径流 量 呈 显著 减少 趋势 。 通 过 以 上 分 析 
结果 可 看 出 ,M-K 检 验 、Spearman 秩 次 相关 检验 结果 
与 线性 趋势 检验 结果 相 吻 合 ,说明 降 水 量 和 径流 量 
趋势 判断 可 靠 度 强 。 尽 管 降水 量 是 径流 量变 化 的 
支配 性 因素 ,但 针对 黄河 上 游 降 水 量 增加 但 径流 量 
显著 减少 这 一 现象 说 明 影响 黄河 上 游 径流 量变 化 
除了 降水 量 这 一 因素 外 ,还 有 人 类 活动 .气候 变化 
等 其 他 因素 。 

3.1.2 黄河 上 游 降水 量 与 径流 量 突 变 分 析 对 图 
3a~b 所 示 的 黄河 上 游 降水 量 和 径流 量 M-K 突 变 检 
验 曲 线 进行 分 析 , 由 i 和 如 曲线 交点 大 致 判断 , 黄 
河上 游 降 水 量 可 能 发 生 的 突变 年 份 为 1962、1976、 
2003 年 ,径流 量 突变 年 份 为 1990 年 。 在 Pettitt 检 验 
下 ,黄河 流域 上 游 降 水 量 在 2003 年 发 生 突变 ,其 变 
异 年 份 所 对 应 的 io=182 ,统计 量 己 .=0.88 ,突变 不 显 
著 ; 径 流量 在 1986 年 发 生 突变 ,其 变异 年 份 所 对 应 
的 ko=7090 ,统计 量 忆 .=0, 突 变 显 著 。 采 用 有 序 聚 
类 对 黄河 上 游 降 水 量 和 径流 量 进 行 突变 分 析 ,降水 
量 和 径流 量 离 差 平方 和 的 时 序 变化 如 图 4 所 示 , 根 
据 图 4a 得 出 降水 量 离 差 平方 和 在 2003 年 达到 低谷 
值 , 因 此 认为 2003 年 为 降水 量 最 佳 变异 分 割 点 ,从 


表 1 降水 量 与 径流 量变 化 趋势 统计 检验 值 


Tab.1 Variation trend of precipitation and runoff inspection values 


类 型 Z 值 Zan 趋势 性 显著 性 r, Im T. 趋势 性 显著 性 
降水 量 0.78 1.96 增加 不 显著 0.09 0.78 2.01 增加 不 显著 
径流 量 -3.82 1.96 减少 显著 -0.51 2.58 2.01 减少 显著 


平 临 界 值 。 下 同 。 


注 :Z 值 为 M-K 检 验 统计 值 ;Zi 为 M-K 检 验 95% 置 信和 度 水 平 临界 值 ;为 Spearman 秩 次 相关 系数 ;| 了 为 1 检验 统计 量 ; 7 为 1 检验 95% 置 信和 度 水 
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图 4 黄河 上 游 降水 量 与 径流 量 离 差 平方 和 


Fig. 4 Sum of squares of the precipitation and runoff in the upstream of the Yellow River 


图 4b 看 出 ,1986 年 为 径流 量 离 差 平 方 和 的 最 小 值 ， 
即 认为 1986 年 为 黄河 上 游 径流 量 突变 年 份 。 采 用 
累积 距 平 对 黄河 上 游 降 水 量 和 径流 量 突变 年 份 进 
行进 一 步 精确 识别 ,根据 图 5 所 示 的 黄河 上 游 降 水 
量 与 径流 量 累 积 距 平 曲线 看 出 ,降水 量 和 年 径流 量 
均 在 研究 时 段 内 发 生 了 一 次 较为 明显 的 突变 ,降水 
量 发 生 在 2003 年 ,径流 量 发 生 在 1986 年 。 综 合 以 
上 分 析 ,确定 黄河 上 游 降水 量 突变 年 份 为 2003 4E , 
1956—2002 年 平均 降水 量 为 415.7 mm ,2004 一 2017 
年 平均 降水 量 为 435.1 mm, 突变 前 后 变化 率 为 
4.67% ;径流 量 突变 年 份 为 1986 年 ,1956 一 1985 年 
平均 径流 量 为 2.49x10" m’,1987 一 2017 年 平均 径流 
量 为 1.61x10" m? ,突变 前 后 变化 率 为 35.34%。 

3.2 分 区 降水 量 与 径流 量 特征 

32.1 分 区 降水 量 、 径 流量 变化 趋 荔 ” 对 各 区 域 搬 
值 的 空间 降水 量 及 径流 量 数据 进行 分 析 , 从 表 2 所 
示 的 分 区 降水 量 和 径流 量 特征 值 可 看 出 ,黄河 上 游 
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图 5 黄河 上 游 降水 量 与 径流 量 累 积 距 平 


Fig.$ Precipitation and runoff accumulative anomalies in 


the upstream of the Yellow River 


年 均 降水 量 由 西向 东 呈 减少 趋势 ,3 个 区 域 之 间 变 
化 率 分 别 为 15.37%、39.53%; 年 均 径 流量 变化 量 由 
西向 东 也 呈 减 少 趋势 , 亚 区 径流 量变 化 量 为 负 值 ， 
说 明 该 区 域 产生 的 径流 量 被 全 部 消耗 。 

根据 图 6a~f 所 示 的 各 区 域 降水 量 和 径流 量变 
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表 2 分 区 降水 量 与 径流 量 特征 值 


Tab.2 Characteristic values of precipitation and runoff in different regions 


类 型 研究 区 域 均值 最 大 值 最 大 值 出 现年 份 最 小 值 最 小 值 出 现年 份 极 值 比 
降水 量 /mm I 区 530.90 658.20 1989 418.60 1956 1.57 
[I 区 449.30 622.50 1967 338.20 1991 1.84 

TIL [X 271.70 425.30 1961 148.30 1965 2.87 

径流 量 /10: m? 工区 198.64 328.60 1989 105.75 2002 3.11 
工区 97.29 214.62 1964 13.67 2009 13.67 

TIL [X -91.80 -58.92 1956 -121.14 2005 0.49 


化 过 程 线 可 知 , I K5 IKKE. mm: 5004 3.5x10? n? (103), 11.5X 10? m°» (103) FI 4.4x 
(10a) 2.8 mm (10a) ' 的 速率 增加 ,而 五 区 以 3.1 10° m°- (102) ARERI o 


mm: (10a)! 的 速率 减少 ; T M 、 焉 区 径流 量变 化 分 采用 M-K 和 Spearman 秩 次 相关 检验 方法 验证 
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图 6 不 同 分 区 降水 量 与 径流 量变 化 特征 


Fig.6 Variation characteristics of precipitation and runoff in different research regions 
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1956 一 2017 年 时 段 内 3 个 区 域 年 降水 量 与 年 径流 量变 
化 趋势 及 在 置信 水 平 o=0.05 下 的 显著 性 ,结果 如 表 3 所 
示 。 两 种 方法 检验 结果 一 致 , [ 区 降水 量 呈 显著 增加 趋 
势 , 卫 区 时 不 显著 增加 趋势 ,了 区 时 不 显著 减少 趋势, 而 
径流 量 均 呈 减 少 趋 势 ,该 结果 与 线性 趋势 检验 一 致 。 

3.2.2 径流 量 对 降水 量变 化 的 响应 分 析 ”原则 上 降 
水 量 是 影响 径流 量 的 主要 支配 性 因素 ,但 气候 变化 


和 强人 类 活动 等 改变 了 影响 径流 量变 化 的 主要 因 
素 。 因 此 针对 黄河 上 游 3 个 不 同 分 区 ,分 析 降 水 量 
与 径流 量 的 变化 关系 ,确定 降水 量 对 径流 量 的 影响 
程度 。 采 用 Kendall、Spearman 和 Pearson 法 做 相关 
分 析 ,判别 时 段 内 降水 量 对 径流 量 的 影响 ,结果 如 
表 4 所 示 。 

3 种 相关 性 分 析 结 果 表明 ,在 年 尺度 上 ,降水 量 


表 3 不 同 分 区 M-K 和 Spearman 检验 结果 


Tab.3 M-K test and Spearman test results in different research regions 


研究 区 域 ”检验 对 象 Z 值 Zo 趋势 性 显著 性 7 m Tu 趋势 性 SHE 
I 区 降水 量 2.45 1.96 增加 显著 0.30 3.26 2.01 增加 显著 

径流 量 -0.97 1.96 减少 不 显著 -0.13 0.89 2.01 减少 不 显著 

工区 降水 量 1.00 1.96 增加 不 显著 0.13 1.14 2.01 增加 不 显著 
径流 量 -3.74 1.96 减少 显著 -0.49 2.50 2.01 减少 显著 

亚 区 降水 量 -0.66 1.96 减少 不 显著 -0.35 1.97 2.01 减少 不 显著 
径流 量 -4.71 1.96 减少 显著 -0.58 2.78 2.01 减少 显著 


与 径流 量 存在 一 定 的 相关 性 ,根据 相关 系数 的 大 小 
可 知 ,3 种 检验 方法 的 检验 结果 均 显示 降水 量 和 径 
流量 的 相关 性 在 区 域 上 呈 工区 > 亚 区 > 开 区 ,说 明 降 
水 量 对 径流 量 的 影响 在 唐 力 辫 水 文 站 以 上 的 河源 
区 最 大 ,其 次 是 下 河沿 到 头 道 抛 水 文 站 之 间 的 耳 
区 ,降水 量 对 径流 量 影响 最 小 的 区 域 是 唐 力 交 到 下 
河沿 之 间 的 区。 因此 以 该 区 域 为 基准 区 ,建立 降 
水 量 -径流 量 双 累积 曲线 (图 7)。 匡 区 和 于 区 作为 
变化 区 ,将 累积 降水 量 代 和 人 基 准 方程 中 ,得 到 累积 
径流 量 ,并 经 过 反复 推算 得 到 径流 量 的 理论 值 , 通 
过 理论 值 和 实测 值 之 间 的 差异 ,确定 降水 量 因素 对 


表 4 不 同 分 区 降水 量 与 径流 量 相 关 性 检验 结果 
Tab.4 Test results of correlation between precipitation 


and runoff in different regions 
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图 7 基准 区 (工区 ) 降 水 量 径流 量 双 累积 曲线 


Fig.7 Double accumulation curve of precipitation and runoff 


in datum zone (zone I ) 


径流 量 的 影响 和 非 降水 量 因素 对 径流 量 的 影响 , 计 
算 结果 如 表 5 所 示 。 当 以 工区 为 基准 区 时 , TE EXE 
水 量 因 素 对 径流 量 的 影响 达到 25.08% , 非 降水 量 因 
素 为 74.92% ; 亚 区 降水 量 因 素 对 径流 量 的 影响 达 
32.1496 , 非 降水 量 因素 为 67.86%。 


表 5 降水 量 与 非 降水 量 因素 对 径流 量 的 影响 程度 


Tab.5 Influence degree of precipitation and non-precipitation factors on runoff 


研究 区 域 Kendall 检 验 。 Spearman 检 验 。 Pearson 检 验 
IK 0.49 0.68 0.71 
工区 0.15 0.21 0.34 
IRES 0.33 0.48 0.43 
研究 区 域 。 平均 降水 量 /mm d a BS 
IK 449.3 9.73 
[X 2711.7 —9.18 


F 均 径流 量 降水 量 因素 对 非 降水 量 因 素 对 
/10" m 径流 量 的 影响 /% 径流 量 的 影响 /% 

1.73 25.08 74.92 

1.05 32.14 67.86 
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4 讨论 


(D 众多 学 者 对 黄河 上 游 径流 量 和 降水 量 的 研 
究 基本 都 集中 在 唐 乃 雍 以 上 的 河源 区 (工区 ) ,对 头 
道 拐 以 上 整个 上 游 地 区 的 研究 较为 缺乏 。 本 文采 
用 黄河 上 游 控 制 站 头 道 抛 水 文 站 降水 量 和 径流 量 
数据 进行 研究 ,但 黄河 上 游 地 处 青藏 高 原 .内 蒙古 
高 原 和 黄土 高 原 交 接地 带 ,控制 面积 较 大 , 仅 采 用 
头 道 拐 控制 站 数据 得 出 的 研究 结论 并 不 能 代表 整 
个 黄河 上 游 地 区 ;同时 年 内 降水 强度 对 人 研究 结果 也 
有 较 大 影响 ,不 同 降水 强度 对 于 径流 量 的 影响 有 往 
进一步 深入 人 研究 。 

(2) 本 文 基于 水 资源 二 级 区 正式 图 件 , 综 合流 
域 自然 环境 和 人 类 活动 对 径流 量 的 影响 ,以 唐 乃 
交 、 下 河沿 、 头 道 拐 水 文 站 断面 为 界 将 黄河 上 游 划 
分 为 3 个 研究 区 域 ,研究 区 域 划 分 主观 性 强 , 较 多 的 
考虑 了 人 类 活动 对 径流 量 的 影响 ,但 气候 等 因素 对 
于 径流 的 影响 也 不 容 忽视 ,如 何 对 研究 区 进行 更 加 
合理 的 划分 仍 需要 考虑 。 

(3) 唐 乃 亥 水 文 站 以 上 工区 ,径流 量 补给 以 降 
水 量 为 主 , 占 源 区 总 径流 量 来 源 的 95.9% 7 ,如 果 单 
纯 考 虑 降水 量 因素 对 径流 量 的 影响 , 即 可 粗略 认为 
该 区 域 径 流量 为 不 受 其 他 因素 影响 的 天 然 径 流 
量 。 因 此 ,以 工区 为 基准 区 模拟 天 然 径流 量 ,得 出 
开 、 焉 区 非 降水 量 因素 对 径流 量 的 影响 很 大 ,这 与 
径流 量 与 降水 量 的 相关 性 检验 结果 一 致 ,与 工区 内 
水 利 工程 林立 ,分布 着 诸如 龙 羊 峡 、 刘 家 峡 等 联合 
调度 的 大 型 水 库 , 下 区 为 黄河 上 游 发 达 的 农业 灌溉 
区 ,人 类 活动 成 为 影响 下、 焉 区 域 径流 量 的 主要 因 
素 这 一 实际 情况 也 较 符 合 ,此 外 也 得 出 了 开 区 非 降 
水 量 因素 对 径流 量 的 影响 大 于 亚 区 的 结果 。 要 更 
加 关注 大 型 水 利 工程 的 修建 与 农业 灌溉 引水 对 径 
流量 的 影响 ,尤其 是 流域 内 大 型 水 利 工程 的 联合 调 
度 对 径流 的 影响 。 尽 管 研究 结论 表明 基于 唐 力 北 
水 文 站 以 上 工区 内 人 类 活动 对 该 区 域 的 影响 非常 
小 ,其 径流 为 无 人 类 活动 干扰 的 天 然 径 流 ,但 此 区 
域 并 非 完 全 没有 人 类 活动 ,例如 ,2011 年 建成 运行 
的 日 调节 水 库 班 多 电站 对 此 区 域 的 径流 产生 的 影 
啊 。 另 外 ,全 球 变 暖 引发 了 冻 土 的 变化 ” ,河源 区 冻 
土 的 退化 也 会 对 径流 产生 影响 ,因此 非 降水 量 因 
素 对 该 区 域 径流 产生 的 定量 影响 仍 需 进一步 解决 。 


(1) 整个 黄河 上 游 ( 头 道 拐 以 上 区 域 )1956 一 
2017 年 降水 量 呈 非 显 著 增 加 趋势 ,倾向 率 为 1.80 
mm (102) ,径流 量 呈 显著 减少 趋势, 倾向 率 为 1.94x 
10* m…(10a)', 二 者 变化 并 不 一 致 。 

(2) 根据 各 类 检验 方法 综合 得 出 黄河 上 游 
1956—2017 年 降水 量 序列 的 突变 年 份 为 2003 年 , 径 
流量 突变 年 份 为 1986 年 。 

(3) 对 黄河 上 游 进行 分 区 研究 , E. TE M 
水 量变 化 分 别 呈 显著 增加 不 显著 增加 和 不 显著 减 
少 趋势 ,而 径流 量 呈 现 非 显著 减少 .显著 减少 和 显 
著 减 少 趋势 。 

(4) 分 别 对 3 个 区 域 的 降水 量 径流 量 进行 相关 
检验 ,结果 表明 ,3 个 区 域 降水 量 和 径流 量 相关 关系 
tg T EX TIE EX TE EX , 唐 乃 净 水 文 站 以 上 的 河源 区 降 
水 量 对 径流 量 的 影响 最 大 ,下 河沿 到 头 道 抛 水 文 站 
之 间 的 陡 区 次 之 , 唐 乃 辫 到 下 河沿 之 间 的 卫 区 降水 
量 对 径流 量 影响 最 小 。 

(5) 以 人 类 活动 影响 较 小 的 工区 为 基准 区 , 通 
过 双 累 积 曲 线 分 析 得 出 工区 降水 量 因 素 对 径流 量 
的 影响 达到 25.08% , 非 降水 量 因素 为 74.92% , M DX 
降水 量 对 径流 量 的 影响 程度 为 32.14% , 非 降 水 量 
素 为 67.86%。 
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Changes of runoff and precipitation in the upstream of Yellow River during the 
past 60 years 


CHENG Yi", WU Lanzhen, LIU Fenggui, SHEN Yanjun’? 
(1. College of Geosciences, Qinghai Normal University, Xining 810008, Qinghai, China; 2. Center for Agricultural Resources 


Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang 050022, Hebei, China; 


3. Water Conservancy and Hydropower Engineering College, Gansu Agricultural University, Lanzhou 730070, Gansu, China) 


Abstract: In this research, linear trend test, M-K trend test, Spearman’ s rank correlation test, M-K test, Pettitt’ s 
nonparametric test, orderly clustering analysis, the cumulative departure method, and the double-mass curve were 
applied to study the precipitation and runoff variation characteristics of upstream of the Yellow River, China. 
Further, the relationship between runoff and precipitation is discussed. Results show an insignificant increasing 
trend in rainfall in the upstream of the Yellow River, while runoff exhibits a distinct decreasing trend. Both 
precipitation and runoff present a significant deflection point, whose change rate between before and after is 
4.67% (2003) and 35.34% (1986), respectively. For the three subregions, subregion I: above Tangnaihai station, 
subregion II: from Tangnaihai to Xiaheyan, and subregion III: from Xiaheyan to Toudaoguai, rainfall presents a 
significant increasing trend, an insignificant increasing trend, and a significant decreasing trend. However, runoff, 
for all subregions, exhibits a significant decreasing trend. With subregion I taken as the baseline, the relative 
impact of precipitation on runoff in subregion II is 25.08%. Consequently, the nonprecipitation factors account 
for 74.92%. Similarly, the percentage for subregion III is 32.14% and 67.8696, respectively. This research will 
offer theoretical references for water resource management and planning. 

Key words: the upstream of Yellow River and its regions; runoff; precipitation; variation characteristics; re- 
sponse relationship 


